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J> . symmetry in the minimizer of the GP functional is broken or not. For 

- —i 

the absolute ground state, we prove that in the GP limit a modified 
%^ ' GP functional depending on density matrices correctly describes the 

energy and reduced density matrices, independent of symmetry break- 
ing. For the bosonic ground state this holds true if and only if the 
symmetry is unbroken. 



We investigate the ground state properties of a gas of interacting 
particles confined in an external potential in three dimensions and 
subject to rotation around an axis of symmetry. We consider the so- 
called Gross-Pitaevskii (GP) limit of a dilute gas. Analyzing both the 
absolute and the bosonic ground state of the system we show, in par- 
ticular, their different behavior for a certain range of parameters. This 
parameter range is determined by the question whether the rotational 



1 Introduction 

Since the first experimental realization of Bose-Einstein condensation (BEC) 
in dilute gases of alkali atoms [U |2] , much interest has been devoted to the 
study of their rotational properties. Beautiful results showing the appearance 
of vortices and the formation of vortex arrays have been obtained in various 
experiments [3J E] . Many of the striking features of rotating Bose-Einstein 
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condensates are well described by means of the Gross-Pitaevskii (GP) func- 
tional, and most theoretical investigations rely on the approximation made 
in its use (see, e.g., jSJE])- 

Recently a rigorous justification of the GP functional for non-rotating 
systems has been obtained |7], and a proof of BEC in the ground state was 
given in [Sj. These results do not extend to the case of a rotating system 
in a simple way, however. In this paper, we investigate the effect of the 
rotation. Several new features come into play. First, Bose statistics become 
essential. While for the case of a non-rotating system the ground state of 
the Hamiltonian is automatically symmetric in the particle coordinates, this 
in not necessarily the case for rotating systems and, in fact, is shown not to 
be the case for the system under consideration here (at least for a certain 
parameter range). That is, it is important to distinguish the absolute ground 
state from the bosonic ground state (the ground state of the Hamiltonian 
restricted to totally symmetric wave functions), the two having significantly 
different physical properties. Secondly, the appearance of vortices breaks 
the rotational symmetry, leading to non-uniqueness of the minimizer of the 
GP functional (see Theorem Q below). This makes it necessary to study a 
generalized GP functional, which depends on density matrices rather than 
densities alone. These two properties are in fact related, as we show, the 
symmetry being broken if and only if the absolute and bosonic ground state 
energy differ by a significant amount (of the order of the energy itself). 

The main result of this paper is that the generalized GP density matrix 
functional mentioned above correctly describes the absolute ground state of 
the Hamiltonian under consideration here, in a certain dilute limit. More 
precisely, it gives the correct asymptotics of the ground state energy and all 
the n-particle reduced density matrices of the ground state of the system. 
In the case of unbroken rotational symmetry, the GP density matrix func- 
tional agrees with the usual GP functional, and these results apply also to 
the bosonic ground state. This extends previous results on superfluidity in 
In the general case of broken symmetry, however, we are not able to 
establish the precise asymptotics in the bosonic case, but we give bounds on 
the corresponding energy that show in particular that it's energy differs from 
the absolute ground state energy by a significant amount. Moreover, this 
implies that the absolute ground state has a huge degeneracy that increases 
exponentially with the particle number. 

Before we can give the precise formulation of the preceding statements, we 
have to define the system under consideration. For f2 6 IR 3 and a realvalued 
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V G L~ (R 3 , dx) let H be the one-particle Hamiltonian 

H = -A-fl-L + V(x) , (1.1) 

acting on L 2 (R 3 , <icc). Here A = V 2 is the Laplacian on R 3 , and L = —ixAV 
denotes the angular momentum operator. The operator Hq is the appropriate 
Hamiltonian for one particle confined in a trap potential V and rotating with 
angular velocity f2, in the rotating frame. To ensure that H is semibounded 
from below, we assume that 

G L°°(R 3 ) (1.2) 

for all uj in some non-zero interval [0, Q c ) C R + , where we denote by [•]_ 
the negative part, and r denotes the radial distance perpendicular to f2, i.e., 
r = \x A fi|/|f2|. We call the largest possible Q c the critical angular velocity, 
allowing it to be infinity. Writing H as 

Ho = HV - A) 2 + V(x) - i^r 2 , (1.3) 

where A = ~f2 A x is the vector potential of a constant magnetic field in the 
direction of f2, we see that Ho is semibounded from below for |f2| < Q c , and 
unbounded for |f2| > Q c . We will henceforth always assume that |f2| < Q c . 
Note that the condition (J1.2)) on V, as well as Q c , depend on the direction 
of f2, which we assume to be fixed throughout the paper, whereas its absolute 
value is allowed to vary. 

We also assume that limia-^oo V(x) = oo uniformly in all directions, 
implying discrete spectrum of Hq. Moreover, we demand that V commutes 
with Q, ■ L, or, in other words, it is axially symmetric, being a function that 
depends only on r and z = x ■ f2/|f2|. Without loss of generality V > 0. 

To define the iV-particle problem, consider the Hilbert space TLn = 
(£)f =1 L 2 (M. 3 ,dXi), the iV-fold tensor product of the one-particle space. The 
Hamiltonian corresponding to particles in a trap V, rotating with angular 
velocity f2, and interacting pairwise with an interaction potential v, is given 
by 

H N ^ a = J2 H ° )+ E ^((Xi-x^/a) . (1.4) 

i=l l<i<j<N 



V(x) - ^r 2 
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Here the superscript (i) means that Hq acts on the z'th part in the tensor 
product. The interaction potential v is assumed to positive, spherically sym- 
metric and of compact support. We do not demand it to be integrable, it is 
allowed to have a hard core, which reduces the domain of definition of H^ ^ a 
to wave functions in that vanish whenever two particles are closer to- 
gether then the size of the hard core. The positive parameter a appearing 
in i?jv,n, a determines the range of the interaction v. We assume that v has 
scattering length 1 (see j7j or JUj for a definition), implying that a~ 2 v{x/a) 
has scattering length a. 

We are interested in the ground state properties of H^ n a for large N 
and small a. In particular, the Gross-Pitaevskii limit a ~ N^ 1 will be in- 
vestigated. Its significance comes from the fact that it ensures that the 
contributions to the energy of all the terms in the Hamiltonian (jl.4|) are of 
the same order as iV —>■ oo. Let E^ M denote the ground state energy of 
HN,n,a, i-e., 

E QM (N, Q, a) = inf spec H N ^,a ■ (1-5) 

Besides N, f2 and a it depends on V and v , of course, but these potentials 
are assumed to be fixed once and for all. From the discussion above it is clear 
that E QM is finite for |fi| < Q c , and that E QM (N, ft, a) = E QM (N, -ft, a). 

A different problem is obtained by restricting ourselves to wave functions 
in Hn that are totally symmetric with respect to exchange of two particle 
coordinates, which corresponds to assuming the particles to be Bosons. We 
denote the corresponding ground state energy by E^ JBe (N, f2, a), i.e., 

E^(N,fl,a) = inf [(^f\H N ^) : G P boae H N , ||*|| 2 = l} , (1.6) 

where Pbose denotes the projection onto totally symmetric functions. The two 
quantities (|1.5|) and (jl.6|) a a priori unrelated, except for the trivial inequality 
E QU (N,n,a) < E^ c (N,fl,a). For ft = 0, i.e., the case of a Schrddinger 
operator, it is well known that E QM (N, 0, a) = E^{N, 0, a), but for O ^ 
this need not necessarily be the case. In fact it will turn out that, at least 
for certain values of the parameters, the two quantities are different. 

Before we can state our main results, we have to introduce some func- 
tional which will turn out to be related to the asymptotic behavior of the 
iV-particle problem defined by H^^^ for large iV and small a. We do this 
in the next section. 
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2 Gross-Pitaevskii Functionals 

We define the Gross-Pitaevskii density matrix (DM) functional to be 

£™[ 7 ] = Tr[tf o7 ] + Amg J p,(x) 2 dx . (2.1) 

Here 7 is a one-particle density matrix, a positive trace-class operator on 
L 2 (R 3 ,dx), and p 7 denotes its density. The corresponding ground state en- 
ergy, the infimum of (J2.1j) under the condition Trpy] = 1, will be denoted by 
£ DM (ft,s),i.e., 

E™(n,g) = inf {^[7] : Tr[ 7 ] = l} . (2.2) 

Note that one could equivalently define the ground state energy under the 
subsidiary condition Tr[ 7 ] = N, or any other constant, but this can be easily 
related to (J2.2)) by a trivial scaling of E DM and g. Again it is clear that 
E BM (Q,g) is finite for |f2| < Q c and g > 0, which we will always assume. 

The analogue of the DM functional for a two-dimensional gas was intro- 
duced in [IT] as a generalization of the standard GP functional, which is 
given by restricting £™ to density matrices of rank one. Equivalently, one 
can write 

£&] = (0 W> + 4tt<? f \<j>{x)\ A dx (2.3) 

for functions G L 2 (M. 3 ,dx), and define the corresponding ground state 
energy as 

£ GP (fi,s)=inf{^] : Uh = l} . (2.4) 

Since one- dimensional projections are legitimate density matrices in ()2.1jl . it 
is clear that E DM ({l,g) < E GP (n,g). 

In the case of two dimensions both the functionals £™ and S^ p g were 
studied in [TJJ. Many of the results translate directly to the three-dimensional 
case, with minor modifications, and we merely state them here in the follow- 
ing two propositions, omitting the proofs. 

PROPOSITION 1 (Minimizer of S™). For each < \Q\ < Q c and 

g > there exists a unique minimizing density matrix for \2.1\) under the 
condition Trpy] = 1, denoted by 7$^. This minimizer also minimizes the 
linearized functional 



7 ~Tr[0ff o + 87raC) T ] 



(2.5) 
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(under the same normalization condition), where p^z denotes the density of 
7njg- Moreover, 7^ has finite rank. Its density p^ is a bounded func- 
tion, with p^^{x) < Psi^/{8iTg), where p^i 1 is the chemical potential in DM 
theory, which is the ground state energy of \2.di) . 

Note that for uniqueness g > is essential, since the ground state of H 
may be degenerate, implying non-uniqueness for g = 0. 

Because of uniqueness, we also know that p^i is axially symmetric, with 
symmetry axis fl. Note that the uniqueness of 7^ is non-trivial. From the 
strict convexity of £™ in p 7 it follows only that p^ is unique, i.e., that 
every possible minimizer has the same density. By examining the ground 
state space of the operator appearing in (|2.5|) one then shows that there can 
be only one density matrix with the property that it minimizes (J2.5)) and has 
Pn^g as ^ s density. 

For S^ F essentially the same results as in Prop. ^ are true, except for 
uniqueness. 

PROPOSITION 2 (Minimizers of ££ p g ). For each < \Q\ < Q c and 

g > there exists a minimizing function for \2. under the condition ||0||2 = 
1 . Any minimizer <p fulfills the GP equation 

# o + 87TsH 2 = /^> , (2-6) 
where p^ g = p^ g (4>) ^ s the chemical potential, given by 

pll = E GP (n,g)+4ng J \cf>(x)\*dx . (2.7) 

One might suspect that the two minimization problems (|2.2j) and (J2.4j) 
are equivalent, in the sense that the minimizer of has rank one, i.e, 
7q^ is a one-dimensional projection onto the minimizer of £^ F , which would 
consequently have to be unique (up to a constant phase factor, of course). 
However, for f2 7^ this is not true, at least not for g large enough. The 
following Theorem ^ is the analogue of Theorem 4 and Corollary 2 in ^T] , 
where the case of a two-dimensional system was considered. Since its proof is 
not a simple generalization of the two-dimensional case, it will be given in the 
Appendix. Due to some technical complications it is more difficult to prove 
this theorem in three dimensions, and it will be convenient not to consider 
the most general external potentials V, but to restrict ourselves to a special 
class with sufficiently nice properties. The general class of V's which our 
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proof applies to is rather difficult to characterize, but sufficient conditions 
are easy to state. We assume them for simplicity, they are general enough to 
allow for a quite large class of V's, but they are by no means necessary for 
Theorem to hold, as the proof in the appendix shows. 

For Theorem^ (and only there) we will assume that the external potential 
V fulfills the bounds 

V(x) < const. (1 + r s + \z\ p ) (2.8) 

and 

V(x) > const. (r s + \z\ p ) - const. (2.9) 

for suitable constants (independent of r and z) and for some exponents 2 < 
s < oo and < p < oo. Note that s > 2 is necessary for condition (jl.2|) to 
hold. 

THEOREM 1 (Non-equivalence of £g* and £™). Assume that V 
satisfies the bounds \2.8\) and \2.ty) . For any < |f2| < Cl c there exists a g\n\ 
such that g > g\a\ implies that no minimizer of £^ g is an eigenfunction of 
the angular momentum f2 • L, and consequently the minimizer is not unique. 
Moreover, for g > g\n\, E GP (Q,g) > E DM (Q, g) , and the minimizer 7^ of 
Sflg has at least rank 2. 

For a given direction of f2, we will denote by H C [0, Q c ) x M + the set of 
parameters where £^ p g and £^ are not equivalent, i.e., 

S = {(\n\,g) : E GP (Q,g) > E DM (Q,g)} . (2.10) 

This is the case if and only if the rank of 7^? is greater or equal to two, and 
therefore, by Prop. [H the ground state of H + 8ngp^ is degenerate. Note 
that Theorem ^ states that S is non-empty, at least for external potentials 
satisfying the bounds (|2.8j) and ()2.9|) . 

In the non-rotating case, i.e., f2 = 0, E GP (0,g) and E DM (0,g) are equal 
for all g > 0. This remains true if |f2| is not too large. In fact one can 
show that there exists an Q g > 0, depending on g, such that E GP (Q,g) = 
E DM (n,g) for |H| < tt g [nj. 

Note that both E GP and E DM are concave functions of their parameters 
and, in particular, continuous. Hence H is an open set, being the complement 
of the zero set of the continuous function E GP (Q, g) — E DM (Q, g), and both 
axis g = and |f2| = are not contained in 5. 
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3 Main Results 

With the preliminaries of the previous section in hand, we can now state our 
main results. We are interested in the ground state energy E^ M (N, f2, a) of 
the iV-particle Hamiltonian (|1.4j) . for large N and small a. In fact it turns out 
that a ~ iV -1 is the case of interest. As explained in the Introduction, a ~ 
iV -1 implies that all terms in the Hamiltonian (jl.4j) yield a contribution to 
the ground state energy of the same order as N — > oo. That is, for some fixed 
g > 0, we will set a = gN^ 1 , or, more generally, we will assume that Na — ► g 
as iV — > oo. Moreover, we will also derive results for the corresponding 
ground states or, more generally, for approximate ground states. We call a 
sequence \Pjv G Hn an approximate ground state if H^jvlh = 1 and 

lim (* N \H N>[l>a V N ) E^ M (N, n, a)' 1 = 1 (3.1) 

N— >oo, Na—*g 

for fixed Q and g. Given such an approximate ground state, we define its 
n-particle reduced density matrix by the kernel 

JV 



T^ ) (x 1 ,...,x n ,y 1 ,...,y n ) = ^- Y] [ ] \ dx 

^ 1 jRS(JV-n) .- LJ - 



ttSSjv " M j=n+l 

x *iv(7r(cci, . . . , aijv))^(7r(2/ 1 , ...,?/„, cc„+i, . . . , as^)) , (3-2) 

where Sn denotes the permutation group and tt(xi, . . . , ccjv) is a permutation 
of the iV variables Xi, 1 < i < N. The * denotes complex conjugation. The 
rff defined in (|3.2|) are trace class operators on 7i ni and the normalization 
is chosen such that Tr[T^] = 1. Our main result on the properties of the 
absolute ground state of H N Si a is the following. 

THEOREM 2 (Ground state asymptotics). For given f2 and g let 7njg 
be the unique minimizer of S^f, with corresponding energy E DM (Q, g) . Let 

denote the n-particle reduced density matrix of an approximate ground 
state of Hn,si,h- Then 

Km ^ u (N,n,gN- 1 ) = E™(n,g) (3.3) 
uniformly in g on compact intervals in (0, oo), and, for each nGN, 

lim rS? = 7 ™ ® ■ ■ ■ ® 7 ™ (3.4) 
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in the usual norm of trace class operators on (££)" =1 L 2 (IR 3 , dxi). 

This theorem will be proved in the next section. Concerning the bosonic 
ground state energy E® oae (N, ft, g), we cannot give the precise asymptotics 
as in (J3.3j) . but we can give upper and lower bounds. One might conjecture 
that (|3.5j) below holds as equality, but we cannot prove this. However, the 
conjecture is supported by the following theorem. 

THEOREM 3 (Asymptotics for bosonic ground state energy). For 

fixed fl and g we have that 

limsupijSgi^n^iV- 1 ) < E GP (n,g) . (3.5) 

Moreover, 

liminf ^ESLiN^gN- 1 ) > E™(n,g) (3.6) 

N— >oo iv 

if and only if 

E GP (n,g)>E™(n,g), (3.7) 

i.e., (\Q\,g) G S, and the same is true with liminf replaced by limsup in 
Eg. Jig). 

Remark 1 . Inequality ()3.6|) holds uniformly in g for compact intervals of g 
in 5. More precisely, for any e > such that the closed interval [(|f2|, g — 
e), (|fi|, g + e)] is contained in H, 

liminf inf l^E^ e (N, ^g'N' 1 ) - E™(n,g')\ > . (3.8) 

This property will be important in the proof of Corollary ^ 

Remark 2 . Theorems 121 and E2 together have the following consequences on 
the Bose gas. If (|f2|,gr) E, e.g., E GP (fl, g) = E DM (Q, g), there is a unique 
minimizer of the GP functional ()2.3|) . Moreover, 

^^(N^gN- 1 ) = E GP (n,g) , (3.9) 

as can be seen from the lower bound ()3.3|) and the upper bound ()3.5|) . There- 
fore the bosonic ground state is an approximate ground state for the unre- 
stricted problem, and hence (J3.4)) holds. The density matrix 7^ in this case 
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is the one-dimensional projection onto the minimizer of£g, and flU proves 
complete Bose-Einstein condensation of all n-particle density matrices. For 
Q = this was proved in jH] . The persistence of BEC for f2 7^ can also be 
interpreted as a superfluid behavior of the system (see 0). 

Note that these results are true not only for the bosonic ground state, but 
for any approximate bosonic ground state. Moreover, the notion of approxi- 
mate ground states is readily generalized to iV-particle density matrices. In 
particular, the assertions above are true for all Gibbs states 

r /3 _ P hose exp(-[3H N ^ a ) 

1 N.fl.a — rp r D I qtt M ' \ 6 ' W 

' ' TT[P hose exp(-/3H N:nA )\ 

with /3 > 0, assuming that the trace is finite, which is guaranteed for external 
potentials V that increase at least logarithmically in \x\ at infinity. We will 
show in Sect. 14.31 that 

for all (3 > 0, implying, for (|fi|, g) ^ 5, complete BEC of the reduced density 
matrices of (|3.10jl . i.e., they converge to the right hand side of ()3.4|) . where 
7^ is now the projection onto the unique GP minimizer. Note, however, 
that fixing f3 is really a zero-temperature limit, since the relevant temperature 
scale depends on the (mean) density, which goes to infinity in our limit. To 
obtain a true effect of the temperature one has to scale it appropriately 
with N. E.g., for a harmonic trap potential the relevant temperature scale 
would be T ~ N 1 ' 3 0. 

Note that for (|3.11|) to hold true it is essential to restrict oneselves to 
the bosonic subspace in (j3.10|) . Without this restriction ()3.11j) (with E^ e 
replaced by will not be true, as can be seen from the non-interacting 

case g = 0. 

A simple corollary of Theorem El is the non-uniqueness of ground states 
of the Hamiltonian H^^ a . Even more is true, namely the ground state 
degeneracy grows exponentially with N. 

COROLLARY 1 (Ground state degeneracy). For (\£l\,g) E E and N 

large enough, no ground state of H^^^^-i has bosonic symmetry. Moreover, 
if j\f(N, \Q\,a) denotes the multiplicity of the ground state of i7jv,n,o; an d 
g) E S, then 

liminf -j- lnJVpV, IfJl^iV" 1 ) > . (3.12) 
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Remark 3. Corollary ^ states that the degeneracy of the ground state of 
HjffagN- 1 grows at least exponentially with N if (|f2|,g) e H. If the external 
potential V does not grow too slowly at infinity, it is not difficult to see 
that in fact j\f(N, gN" 1 ) grows also at most exponentially. Namely, if 
we assume that Tr[exp(— /3H )] < oo for some (3 > 0, which is guaranteed if 
V has at least logarithmic increase in |a;|, then 

limsup — lnJVXW, \nigN~ 1 ) < lnTr[exp (f3{E DM {tl, g) - H ))] . (3.13) 

This is shown in Sect. 14.41 

The proofs of the results stated in this section can be found in Section 0] 
below. The energy difference of the bosonic system to the one without sym- 
metry restrictions can be understood as being due to the positive correlation 
energy that gets added when trying to symmetrize a state where not essen- 
tially all the particles occupy the same state, as is the case in the regime of 
symmetry breaking, where the minimizer of £ DM has at least rank 2. It may 
therefore be favorable for all bosons to occupy the same state, restoring the 
complete BEC but breaking the rotational symmetry by choosing one of the 
minimizers of the GP functional. 

All the results in this paper refer to three-dimensional systems, but anal- 
ogous results can be obtained also for two-dimensional systems, as has been 
shown, in the non- rotating case, in jT5j. The DM and GP functionals in this 
case have been studied in [TT] . Moreover, we could also allow for internal 
degrees of freedom of the particles, e.g. spin, which would not affect the 
results on the absolute ground state, but the ones on the bosonic ground 
state. In this case the GP functional has to be replaced by the DM func- 
tional restricted to density matrices of a definite rank equal to the number 
of internal states (compare with Section 6 in [TT]). 

4 Proofs 

Before we give the proof of the results stated in the previous section, we state 
in Section 14.11 two auxiliary Lemmas that will be used later. The proof of 
Theorem |21 is given in Section |4~21 and the proof of Theorem [21 together with 
the assertions made in Remarks 1 and 2, in Section l4~3l Finally, Corollary Q 
and the subsequent Remark 3 are proved in Section 14.41 
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4.1 Preliminaries 



The following Lemma is needed in the proof of Theorem|2]below. It is related 
to the generalized Poincare inequalities studied in For a measurable set 
O C M 3 we denote by O c its complement, and by \0\ its Lebesgue measure. 



LEMMA 1. Let A e L^ oc (IR 3 ; 
lirni a .i_ >0O V(x) = oo. Let 

E = inf spec (-zV - A(cc)) 2 + V(x) 



and V e L^ C (M 3 ;1R) ; and assume that 



(4.1) 



and let P denote the projector in TC 
ground states. Let 



dx) onto the corresponding 



AE = inf spec [(-iV - A(x)f + V{x)] t(i-p)w - ^ 



(4.2) 



denote the gap in the spectrum above the ground state energy, which is positive 
because of the discrete spectrum of the operator under consideration. 

For all e > there exists a 5 > such that for all Oct 3 with \O c \ < 5 
and for all f EH 



(, V + A)/| 2 + / |(*V + A)/| 2 + 



v\f\ 2 

>E\\f\\l Hm +AE\\f-Pff LHm . (4.3) 



Proof. It is no restriction to assume that V > 0. As in [13J we will use 
a compactness argument. Suppose that the Lemma is wrong. Then there 
exists an e > and a sequence of pairs (f n , O n ), such that lim^oo \ 0^\ = 0, 
H/nlU^wsi = 1, and 



lim 



e [ |(zV + A)/ n | 2 + / |(*V + A)/ n | 2 

+ / V\f n \ 2 -AE\\f n -Pf n \ 



< E . (4.4) 



so we can 



Now both (iV + A)f n and f n are bounded sequences in L 2 (I 
pass to a subsequence that converges weakly in L 2 (M 3 ) to (iV + A)f and 
/, respectively. We may also assume that \0^\ is finite. Defining T, N 
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by £ 



N 



\ Un>Af ®n we have Sat C O n for n > N. Using weak lower 



semi continuity of the norms in question, we therefore get 



lim inf 

n— >oo 



zV + A)/ n | 2 + / |(zV + A)/ n | 2 + / V\f, 



> sup 

N 



(W + A)f\ 2 + I |(zV + A)/| 2 



[l+e ) f \(^ + A)f\ 2 +[ V\f\ 2 }>E[ |/| 2 . (4.5) 

Now since V goes to infinity at infinity, —A + V has a compact resolvent 
(cf., e.g., [TIJ Thm. XIII. 65]), and hence also (iV + A) 2 + V has a compact 
resolvent [T51 Thm. 2.7]. Since 



/ (\( t V + A)f n \ 2 + V\f n \ 2 )<C 

JWL 3 



(4.6) 



for some C < oo independent of n, we can conclude that f n is contained in a 
compact subset of L 2 (R 3 ), and thus f n —*f strongly in L 2 (R 3 ). This implies 
that \\fW2 = 1, and also 



(4.7) 



lim \\fn-Pfn\\L^) = Wf-Pfh* 
n— >oo 

Together with ()4.5|) this contradicts (|4.4|) . 



□ 



Using more sophisticated methods, as in ^3J, it is possible to investigate 
the relation between e and 5. This is needed to get precise error estimates, 
but we shall not do this here. 

In the proof of Theorem |21 it will be necessary to study a homogeneous 
gas of n particles, described by the Hamiltonian 



n 1 



(4.8) 



l<j<j'<n 



acting on L 2 (A n ). The particles are confined to A = [0,L] 3 , a box of side 
length L, and we use Neumann boundary conditions. A lower bound to the 
ground state energy of (|4.8j) was obtained in jXB], and we shall describe the 
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result here. In fact, the lower bound in ^B] was obtained via a lower bound 
to the expression 

Ef|/ |VJ| 2 + (1-.) / |Vjf]+W - 2 v(( Xl - Xj )/a)\f\\ 

(4.9) 

where < e < 1 and A'j C A™ is given by 

Ki = < (x 1} . . . , x n ) E A n : min \x { - x k \ < R > (4.10) 

for some R > 0. This is exactly the expression that we have to bound from 
below in the proof of Theorem El see Eq. (|4.2fij) . The result is the following. 
It is valid for all spherically symmetric v > with finite range and scattering 
length 1. 

LEMMA 2. Let E e (n,L) denote the infimum of ^4-9j) over all functions 
f e L 2 (A n ) with ||/|| 2 = 1, and let Y = a 3 n/L 3 . If R > aY~^ 17 , e > Y 1 ' 17 , 
and n > Y~ x l 17 , then there exists a constant C > such that 

E £ (n,L) > Arca^ (l - CY 1 ' 17 ) (1 - e) . (4.11) 

The proof can be found in ^Bj (see also [Tj\ for a more elaborate dis- 
cussion). Strictly speaking, it was derived for the expression (|4.9|) with e/2 
replaced by e in front of the first term, but it is easy to see that this addi- 
tional factor does not affect the main result, only the constant appearing in 
(|4.11|) . Note that the condition on n means that n > (L/a) 1 ^ 6 . 

4.2 Proof of Theorem [2] 

We start by considering the case n = 2 in (|3.4|) . In order to be able to obtain 
information about the two-particle density matrix, we start by introducing a 
modified DM functional. Let U G C^°(IR 3 x M 3 ) be realvalued and symmetric, 
i.e., U (x, y) = U (y, x), and let 8 £ R. The modified DM functional is defined 

as 

S DM [ 1 } = Ti[H 0l } + 47cg j Pj {x) 2 dx + 5 J U(x,y)p 1 (x)p 1 (y)dxdy . (4.12) 
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We suppress the dependence on the parameters for simplicity of notation. 
We assume that 5 is small enough such that £ DM is still strictly convex in 
p 7 , which is in particular the case for <5||£/||l2( R 6) < Aug. The ground state 
energy of (|4.12|) will be denoted by E DM , which depends, for fixed U, on f2, 
g and 5. As for the DM functional with 5 = 0, one can use standard methods 
to show the existence of a minimizing density matrix for £ DM . Uniqueness is 
not clear, however, but this is of no concern to us. (For 5 small enough, £ DM 
is strictly convex in p 7 , and therefore the density of a minimizer is unique. 
This convexity property will be important in the proof of Lemma 0] below.) 
Let 7 denote a minimizer of £ DM , with corresponding density p. One can 
show that p is a bounded, continuously differentiable function that decreases 
exponentially at infinity. Moreover, as in Prop.^ 7 also minimizes the linear 
functional 

7 ^ Tr[# 7 ] , (4.13) 

with H given by 

H = H + Sngp(x) + 25 J U(x, y)p(y)dy . (4.14) 

That is, the range of 7 is contained in the span of the ground states of H, 
whose ground state energy is given by 

inf spec H = E BM (n,g,5) +Ang J p(x) 2 dx + 6 J U(x,y)p(x)p(y)dxdy . 

(4.15) 

It is important to note that E DM is differentiable in 5 at 5 = 0. This 
follows from concavity in 5 and the fact that the minimizer at 5 = is 
unique (cf., e.g., [IE])- The derivative is given by 



dE DM (Q,g,6) 



05 



5=0 



J U(x,y)p™(x)p™(y)dxdy . (4.16) 



We will now evaluate upper and lower bounds on the ground state energy 
of the iV-particle problem. We also add our auxiliary potential U to the 
Hamiltonian, and estimate in the following the quantity 



E QM {N,n,a,5) = inf spec 



25 



^ „x 5i 

l<i<j<N 



(4.17) 
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We start with the upper bound. It is derived by analogous considerations 
as in [2j , but instead of using a trial vector we use an iV-particle trial density 
matrix, whose kernel has the form 



N 



T(x 1 , ...,x N ,y 1 ,...,y N ) = F(x u . . . , x N )F(y 1 , , y N ) JJ *y(x u y { 



i=l 



(4.18) 

where 7 is a minimizer of ([4. 12)1 . and F is the Dyson wave function defined 
in [Jj. It is given by 



N 



i=i 



where £3 = min{|£Cj — Xj\, 1 < j < i — 1} is the distance of X{ to its nearest 
neighbor among the points Xi, . . . , JCj_i, and / is a function of t > 0. It is 
chosen to be 

m _{Mwm ^ : <i ( , 20) 

where fo is the solution of the zero energy scattering equation for the interac- 
tion potential, i.e., — A/ + \a~ 2 v(x/a)f = 0, and b is some cut-off parameter 
of order b ~ N" 1 ^ 3 . The function F is a suitable generalization of the func- 
tion Dyson used in ^H] to obtain an upper bound on the ground state energy 
of a homogeneous Bose gas of hard spheres. The calculation of the upper 
bound follows along the same lines as in [7j, with the result that 

E QM (N,fl,gN-\5) < NE DM (n,g,5)(l + 0(N~ 2/3 )) , (4.21) 

uniformly in g on compact intervals. 

We now proceed with the lower bound. Fix some R > and < e < 1. 
They will be chosen later to depend on N in a definite way. We introduce 
the short hand notation 

Xi = (xi, . . . , xt-i, x i+1 , . . . , x N ) (4.22) 

and 

N 

dXi = Yl dx i ■ ( 4 - 23 ) 
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Moreover, let A(x) = ~fl x x, and let rj denote the radial distance of Xi 
orthogonal to f2. Let Vj = V^.. Given any \1/ G T~Cn with ||^||2 = 1, we can 
write the expectation value of our modified Hamiltonian as 



25 



i<j 



\ N 



(4.24) 



where 



*3 



(i) 



(JV-l) 



(1-e) / dcc i |(iVy + A( Xj ))y\ 
+ ~ I ^1(1^ + Ate))*! 2 

1 J9? 



d Xj [ v(xj) - i^Lr| + 87r^pte) + 26 J u ( x v y)p(y) d y ] j M 2 



(4.25) 



and 



E 



(2) 



1-e) [ dx j \(iV j + Ate))^| 2 

+ ~ I dxjliiVj + Aixj))^ 2 

- J ^dxj \ 87rgp(xj)+2S J U(x jiy )p(y)dy^ |*| 2 

N f ( 1 25 \ 

+| V / da;,- f — v(te -Xj)/a) + — U{x i ,x j ) ) |*| 



We choose, for fixed X 



(4.26) 



3 ' 



0,- = < £Cy G M : min \xj — x k \ < R 

k, k=ij 



(4.27) 



In the following, we will investigate the two terms ()4.25|) and ()4.26|) sepa- 
rately. The results are formulated in the following Lemmas. It is always 
understood that, for fixed g, Na — g = o(l) as N — > oo, and that 5 is small 
enough, as explained in the beginning of Section 14.21 
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LEMMA 3. Let P denote the projector onto the ground states of H, defined 
in If R < 0(iY~ 1/3 ) one can choose 0(N~ 2 / 17 ) < e < o(l) as 



N — > oo suc/i t/iat 

AT 



3=1 ^ J 



+5 j U(x,y)p(x)p(y)dxdy + CTr[T^(l - P)]\ (1 - o(l)) (4.28) 

/or some constant C > (depending on CI, g and 5). Here rjp denotes the 
one-particle reduced density matrix of ^ . 

LEMMA 4. // i? > 0(iV- 7 / 17 ) and o(l) > e > 0{N~ 2 / 17 ) as N -> oo ; 
t/ien 

> -iV (4th? y p(z) 2 ^ + 5 ^ U(x,y)p(x)p(y)dxdy^ (l + o(l)) . 

(4.29) 

Before proving these two Lemmas, let us show that they lead to The- 



orem |2j Inserting the lower bounds to i?j and ^ • £jj in (|4.24 
obtain, for an appropriate choice of R and e, 



we 



N 



26 



i<j 



> \E UM (n,g,6) + CTr[T^\l-P)] (1 - o(l)) (4.30) 



as iV —>■ oo, if iVa — ► (7. Together with the upper bound (|4.21|) this implies 
that 

hm ±-E<* M (N,n,gN-\6) = E BU (Sl,g,5) (4.31) 

N— >oo iV 

for all values of Q, g and 5 small, uniformly in g on compact intervals, and 
also that 

lim Tr\T$(l-P)] = . (4.32) 

Eq. 1)4.31)1 for 5 = proves ()3.3j) . Moreover, since E DM is differentiate 111 
5 at 5 = 0, with derivative given in (|4.16j) . we also infer that (see jTHj for 
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details) 

lim f U(x,y)p ( *\x,y)dxdy = [ U(x,y)p™(x)p™(y)dxdy , (4.33) 

(2) (2) 

where p N is the density of the reduced two-particle density matrix T ( N J 
of an approximate ground state at 5 = 0. This is true for all symmet- 
ric U G C^°(IR 3 x IR 3 ), i.e., p$(x,y), being itself symmetric, converges to 
Png{ x )Pn^(v) i n ^ ne sense °f distributions. 

To show the convergence ()3.4j) of the density matrices, not only their 
densities, we proceed as follows. From ()4.32j) at 5 = we infer that 

lim Tr[rSJ } P DM <g> P DM ] = 1 , (4.34) 

N^oo 

where P DM is the projector onto the ground states of H with 5 = 0, which is 
the operator appearing in the linear functional (|2.5|) . Note that P DM is finite 
dimensional, and pDM^MpDM = ^dm by p rop [JJ From (ICTD we infer that 

(2) 

there exists a subsequence of T N that converges to some T in trace class 
norm, with P DM ® P DM r? DM ® P DM = T. 

We want to show that the convergence (|4.33|) necessarily implies that 
T = 7njg ® 7n|g • To do this, we have to take a closer look at the ground 
states of H for 5 = 0. It is clear that they can be taken to be eigenfunctions 
of Q L, and that there is at most one ground state for any given eigenvalue of 
f2 ■ L. For eigenvalue m G Z, the ground state wave function can be written, 
in cylindrical coordinates x = (r, (p, z), as 

f(r,z)e im * , (4.35) 

with f(r, z) > for all r > and all z, and f(r, z)r~' m ' bounded (see JT] for 
details). Note that certainly m > for any ground state. Expanding V in 
terms of these eigenfunctions, and using the fact that because of (|4.33j) its 
density is necessarily given by Pn^^p^^x'), we obtain the equation 

PZ( X )P™ 9 ( X ') = E ^/y(r^)/ fc (r'^ / )^(r^)/ m (r , ^ , )e i(i -° ¥, e i ( fe -" 1 ^' 

j,k,l,m 

(4.36) 

for some coefficients ajkim that determine T. This sum is finite, and since 
the left side does not depend on (p and (p f , we can infer that only terms with 
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j = I and k = m contribute, i.e., ajkim = for j ^ I or k ^ m. Here we 
have also used that terms with the same j — I, but different j and I can not 
cancel, because of the different asymptotics as r — > 0, namely r J+/ . By the 
same reasoning, we have 

^>) = E A ^(^)i 2 ( 4 - 37 ) 

for some < Xj < 1. Using again the knowledge of the r — ► asymptotics 
of the /j, we infer that cijkjk = XjXk, and therefore T = 7^ ® 7n!g- This 
proves ()3.4|) for n = 2. 

It remains to prove Lemmas El and EJ 

Proof of Lemma For fixed Xj define fj by 

fjixj) = ^{x u ...,Xj,...,x N ) , (4.38) 

and let 

W{x) = V(x) - ^r 2 + 8ngp(x) + 25 J U{x, y)p(y)dy . (4.39) 



We have 



with 



Ef >(l- £ ) / dX i F i --^-||[W1. 



(4.40) 



F i= / |(zV + A)/,| 2 + f / |(zV + A)/,| 2 + / (4.41) 
The result now follows by applying Lemma ^ to F,-, noting that \O c A < 



N^-R 3 = o(l) and that 



JV 



£ / dXM - P/J 2 = NTr[T^\l - P)] (4.42) 

jR3(iV-l) 

(compare with 1)4.14)1 and (|4.15)l ). The constant C is then given by the 
spectral gap of H above its ground state energy. Note that the last term in 
(14.40)1 is finite, by our assumption p. 2)1 on V and the fact that U is bounded 
andpeL^R 3 ). □ 
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Proof of Lemma^4\ To obtain a lower bound, we can use the box method, as 
in PE1IZI- More precisely, we divide M. 3 into boxes of side length L, labeled 
by a, and distribute our iV particles over these boxes. Taking Neumann 
boundary conditions in each box and minimizing the energy with respect to 
all distributions of the particles, this can only lower the energy. 
Let 



W a = sup 



4%gp(x)+5 U(x,y)p(y)dy 

Jr 3 



(4.43) 



and 



U aP = inf U(x,y) . (4.44) 

Using the diamagnetic inequality in the two terms in ()4.26|) containing A, 
we can set A = for a lower bound. With E e (n, L) defined as in Lemma El 
we thus obtain 



N 



J2 E , 

i=i 



(2) 



> 



inf \ \^ (E £ (n a , L) - 2W a n a ) + ^7] U a pn a n p - -Y] U aa n a \ , 

{n a } iV iV 

(. a a, /3 a ) 

(4.45) 

where the infimum is taken over all distribution of the N particles into the 
boxes a. Each box contains n a particles, and Y2 a na = The last term 
in (|4.45|) is due to the fact that there are only \n{n — 1) pairs of particles 
in the same box. It can easily be estimated by |5||{7||oo, which is negligible 
compared to the other terms of order N. 

Let n a be a minimizing configuration of the n Q 's in (|4.45|) . Let Y = 
a 3 N/L 3 , and let A,, denote the collection of those boxes a where n a > rjNL 3 . 
If we choose e > Y 1 / 11 and R > a(a 3 i]N)~ 5 ^ 17 (compare with the conditions 
stated in Lemma EJ) we can use ()4.11j) to estimate, for a G A^, 

E E (n a ,L) >4vra^|(l-Cr 1 / 17 )(l-£) , (4.46) 
where we estimated n a by in the error term. For a ^ A^ we simply use 
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E e (n a ,L) > 0. This gives 

(1133) 

'NL 3K JK N 



aeA v ,/3eAr, 



qGA, ; 

WU- (4.47) 



Here e = 1 — (1 — e)Na/g = o(l) for large iV. 

In the following, we will choose L = o(l) and 77 = o(l) as iV — > oo. The 
term in brackets in the last line in (|4.47j) can then be bounded by const. L~ 3 , 
which is, when multiplied by 7]NL 3 , of lower order than N and therefore 
negligible compared to the main terms of order N. 

Consider now the first two terms in (|4.47jl . If 5 is small enough (indepen- 
dent of iV and L) , the term in square brackets defines a positive matrix B a p 
for iV large enough. Moreover, denoting p a = sup xea p(x), 

W a = NL 3 J2 B ^Pp + L , (4.48) 

P 

where £ a = + £a \ with 

4 1 ) = Ang {CY 1 ' 17 + e - eCY 1 ' 17 ) p a (4.49) 

and 

= 5 ^sup J U{x, y)p{y)dy -^11 u ^P^j ■ ( 45 °) 

In the following we denote by B the matrix with coefficients B a p, and likewise 
by W, p, £ the vectors with components W a , p a and £ a , respectively. Let 
also n denote the vector with entries n a for a E A v and for a ^ A^. With 
these notations we can rewrite the first two terms on the right side of (14.47)1 

as 

(n\B\n) - 2(W\n) 

= (n - NL 3 p\B\n - NL 3 p) - (NL 3 p\B\NL 3 p) - 2(£\n) . (4.51) 
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Because of positivity of B, this expression is bounded below by 

-(NL 3 p\B\NL 3 p) - 2N\\£\\ 00 . (4.52) 
Now, if we choose L — o(l) as N — ■> oo, we have 

(iVL 3 ) 2 (p|B|p) = N (attq J p{xf + 5 j U(x, y)p(x)p(y)) (1 + o(l)) 

(4.53) 

because the sums are Riemann sums for the corresponding integrals. More- 
over, ||.£||oo < which is obvious for (|4.49jl . since p Q is bounded, and also 
for (j4.50J) by the same Riemann sum argument as above. This proves the 
Lemma. □ 

This finishes the proof of ()3.4|) for n = 2. The general case n > 1 follows 
in the same manner, perturbing the Hamiltonian (jl.4|) by 

^A^T E U(x h ,...,x in ) , (4.54) 

l<«i<---<i„<V 

where [/ G C(j , °(lR 3n ) is a realvalued symmetric function of n variables, and 
also perturbing the DM functional by a term 

5 U(xi, ...,x n )p y (xx)---Py(x n ) . (4.55) 

JR 3n 

We will only sketch the proof here, pointing out the differences to the case 
n = 2. Following the same steps as in the proof of Lemma 0] above, one sees 
that minimizing over particle numbers in different boxes is, after passing to 
the limit iV — > oo, effectively equivalent to minimizing a functional 



/ [47igp(x) 2 - W(x)p(x)] + 5 / U(xi,...,x n )p(x 1 )---p(x n ) 

(4.56) 



P 

over functions p G L 1 (M 3 ) with p > and j p = 1. Here 

W(x) =8%gp(x) +n5 U(x, x 2 , . . . , x n )p(x 2 ) ■ ■ -p{x n ) , (4.57) 



(n-l) 



and p is the density of a minimizer of the modified DM functional with 
additional term (J4.55)) . with j p = 1. To obtain the analogue of Lemma 0] we 
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have to to show that the minimizer of (J4.56)) is given by p, at least for 5 small 
enough (compare with 1)4.51)1 . which is just a discretized version of 1)4.56)1 for 
n = 2). Compared to the case n = 2 there is an additional difficulty here 
since the functional ()4.56)1 is not necessarily convex. (The convexity in the 
case n = 2 was used in the form of positivity of B.) 

To show that, for small enough 5, the minimizer of the functional (J4.56)) 
is given by p, we rewrite (14.56)1 as 

-Aug J p{xf - {n - 1)5 J U{x u x n )p{x x ) ■ ■ ■ p{x n ) + F[p] , (4.58) 

with 



T\p\ = 4tt 5 J{ P -p) 2 

n „ 

+ y^d (x h ..., Xi) (p(aii) - p(xx)) ■ ■ ■ (p(xi 



p( Xi )), (4.59) 



where 



U {%) (x h ...,Xi)= U(x 1 ,..., x n )p(x i+ i) ■ ■ ■ p(x n )dx i+1 ■ ■ ■ dx n 

(4.60) 

and q 6 N are appropriate integer coefficients. Now T\f\ = 0, so any 
minimizer p must fulfill T[p] < 0. Since ||p — p||i < 2 by assumption, we 
have F[p] > 4:Trg\\p — — const. 5||Z7||c», and hence 

\\p — p\\t ^ const. S (4.61) 
for a minimizer p. On the other hand, 

Hp] > 47r<?||p - p\\l - const. s\\u\\ 2 (\\ P - + ||p - m) , (4-62) 

which is, because of ()4.61)) . positive for 5 small enough, and zero only for 
p = p. Therefore p is the unique minimizer of ([4.56)1 . 

The convergence of the energies for a small interval of 5's around zero 
implies again the convergence of the n-particle densities at 5 = 0, and by 
analogous arguments as in the case n = 2 one can show the convergence of 
the density matrices. The details are left to the reader. 

This finishes the proof of Theorem |21 
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4.3 Proof of Theorem [SI 

We start by computing an upper bound to E^ c . Note that we can not use 
(14. 18(1 as a trial state, for two reasons. First of all, F is not a symmetric 
function and, secondly, the tensor product of one-particle density matrices 
is not a bosonic density matrix, if their rank is bigger than one. The first 
problem turns out not be serious (see below), but the second is. 
Instead of ([4.18)1 . we use as a trial function 



*(aci, . . .,x N ) = F(x u 



x N 



N 

n 

i=l 



sGP, 



(4.63) 



where GP is a minimizer of the GP functional, and F is given by the Dyson 
wave function ([4.19)1 . Note that F is not a symmetric function, but we claim 
that ([4.63)1 nevertheless gives an upper bound to E^ e . More precisely, we 
will show that any trial function of the form ([4.63)1 with F realvalued but not 
necessarily symmetric gives an upper bound to E^ e . To see this, we write, 
for \l/ given in ([4. 63)1 . 



(* \H N ,ci, a \ *> = ^ w> + e / n \0 Gp M\ 2dx 

+ (| f ^v{{ Xl - X ,)/a)-Mg\^{x 3 )A\F 

V i=l, / 



|ViF| 



ni0 Gp (^)i 2 ^ fc , 

k=l 

(4.64) 

where we used partial integration, the GP equation (|2.6|) . and the fact that 
F is real by assumption. The reality of F has the effect the both the cross- 
term in the partial integration and the term with L acting on F vanish. The 
infimum of ()4.64|) over all F is attained for a positive F, and therefore we 
can proceed as in (201 p. 15] to show that the infimum of (|4.64|) over all F is 
the same as the infimum over symmetric F. 

It is therefore legitimate to use ([4.63)1 as a trial function for the bosonic 
ground state problem. The calculation of the expectation value of Hn,si,o, 
follows again along the same lines as in jjj, with the result that 



EfZiN^gN- 1 ) <NE GP (n,g)(l + 0(N-^)) . 



(4.65) 



Letting N — > oo this proves ()3.5[) . 
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From (J3.5)) it follows immediately that 

limsup-i^gi^n^JV- 1 ) > E™(Q,g) (4.66) 

implies (IH.7|) . Since (J3.6)) clearly implies ()4.66j) . it remains to show that 
(COai follows from (j3~7|l . To see this, note that E GP (n,g) > E DM (Q,g) 
implies that the unique minimizer 7™ of has at least rank 2. Therefore 
7ng<8>7ng is noi a bosonic two-particle density matrix, it has both symmetric 
and antisymmetric eigenf unctions. Since the space of bosonic two-particle 
density matrices is closed under the trace norm, ()3.4|) cannot be true for the 
bosonic ground state in this case. Theorem |21 implies that the bosonic ground 
state is not an approximate ground state for the unrestricted problem, and 
hence f!3.6|) holds. This proves Theorem EJ The uniformity of the limit (|3.4j) 
(see the definition of approximate ground state in (j3.1j) ) implies uniformity 
in ()3.6j) . as stated in ()3.8|) in Remark 1. 

We conclude this section with a proof of the assertion (j3.11|) made in 
Remark 2 that the Gibbs states (|3.1U|) are approximate ground states. Let 
T = 1/(3. We start by estimating the free energy 

F(T, N) = -- lnTr[P bosc exp(-(3H N )} . (4.67) 

For simplicity of notation we suppress the dependence on f2 and a from now 
on. It is always understood that a = gN^ 1 . The function F is monotone 
decreasing and concave in T. This implies immediately that F(T, N) < 
F(0,N) = E^ e (N). To obtain a lower bound, we proceed as follows. Let 
n(e) denote the number of eigenvalues of H N \p hosc n N that are smaller or 
equal to e. Then 

(4.68) 

Since v > 0, n(e) is bounded above by the number of eigenvalues of the 
non- interacting Hamiltonian X^-^o tp b0 s C Hjv are sman er or equal to e. 
This, in turn, can be bounded above by 

Tr[Pbo Se exp (ae - «X)i#o )] = exp(ae)C Qi7V (4.69) 



F(T, N) 



(3 



■QM 



(AT) 



n(e) exp(— (3e)de 
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for any a > 0. Choosing a < (3 we therefore obtain 

F(T, N) > ^B^(JV) - i hi g%£) . (4.70) 

Now let e«, i > 0, denote the ordered sequence of eigenvalues of H , including 
degenerate eigenvalues, and let j denote the multiplicity of the lowest eigen- 
value e . The expression C a ^ can easily be bounded above using the fact 
that any totally symmetric eigenfunction of £\ Hq can be characterized by 
the occupation numbers rij of the eigenfunctions of H corresponding to the 
eigenvalue ej. Thus C a> ^ can be expressed as 

C a , N = e - aNeo Tr [p bosc exp(-aEti(^d 4) - e o)) 

e -«E%-fe-eo) ^ (4-71) 



g-aVeo 



where we introduced the notation No = N U {0}. To bound this expression 
from above, we simply neglect the condition Y2j>o n j = N. More precisely, 
we estimate 

N N 

C a>N < e - aNeo J2--- Yl 2 e~ a ^ n ^- e ^ 

™o=0 Wj -i=0 {%£N }j>j 
oo 1 

= (^+l)^-^ n ____. (472) 

Note that the last product is finite, and independent of N. Inserting the 
bound ()4.72|) in ()4.7U|) . letting N oo and afterwords a — > 0, we obtain 

Jim I (F(T, iV) - E^(N)) = . (4.73) 
To see that the same convergence holds for the energies 

E(T, N) = nP^eH N eM-^)} 
Ti[P hose exp(-(3H N )\ 

note that 

dF f T N) = i (F(T, N) - E(T, N)) . (4.75) 

The concavity of F(T, N) in T and (ETTBl imply that N~ l dF/dT -4 as 
iV — > 00, which proves the desired result. 
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4.4 Proof of Corollary Q] 

The first statement on the impossibility of having a totally symmetric ground 
state follows immediately from Theorem El To estimate the ground state de- 
generacy, we proceed as follows. Since the Hamiltonian under consideration 
is symmetric with respect to exchange of two variables, any ground state be- 
longs to a finite dimensional space of an irreducible representation of Sn, the 
permutation group of iV variables. These representations are conveniently 
labeled by Young tableaux. Given the ground state space of Hn^^n- 1 and 
one of the corresponding representations of Sn, let, for fixed (|f2|,g) G H, 
Rn and Cn denote the length of the longest row and column of the Young 
tableaux, respectively. We will show that Cn is bounded above independent 
of N, and Rn/N < 1 — £ for some e > and N large enough. This will allow 
us to show that the dimension of the representation increases exponentially 
with N. 

We first compute an upper bound on Cat. There exists a ground state of 
Hn&^n- 1 which is antisymmetric in the first Cn variables. Neglecting the 
interaction with the other particles, we therefore have 

E QM (N, ft, gN- 1 ) > Tt[H P Cn ] + E QM (N - C N , ft, gN' 1 ) , (4.76) 

where Pc is the projection of H onto the first C eigenvalues. By analogous 
considerations as in the upper bound to E® M (cf. also [3 Thm. III. 2]) it is 
not difficult to show that 

E QM {N, ft, gN' 1 ) - E QM {N - C, fl, gN' 1 ) < const. C (4.77) 

for fixed g and f2, independent of N. Since 

lim ^Tr[H P c ] = oo (4.78) 

because of our assumption that V goes to infinity at infinity, this implies 
that Cn has to stay bounded as iV — > oo. 

We now derive an upper bound on R N . Since there is a ground state of 
^jvn,ffjv-! which is symmetric in the first Rn variables, we have, by analogous 
considerations as above, 

E QM (N, n, gN' 1 ) > E§™(R N , n, gN" 1 ) + (N - R N )E , (4.79) 
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where E = inf spec H . Let rj = limsup^^oo Rn/N, and choose a sub- 
sequence such that Rn/N — > r/ as iV ^ oo. It follows from (|4.79jl and 
Theorem |21 that 



-i 



rj < (E DM (n, g) - Eo) fliminf -£-Eg*(R N , ft, gR~ N l Rn/N)) - E, 

(4.80) 

Since certainly Rn — > oo with N, we can use (13. 8j) to conclude that 

liminf J-^ e (^, ft, gi?- 1 ^/iV)) > £ DM (ft, grj) (4.81) 
if 1 — 77 is small enough such that (|ft|, grj) G H. This shows that 77 < 1 for 

(I^U)eS. 

We have thus shown that, for some e > 0, i?jy < iV(l — e) and Cn < 
1/e for iV large enough. Therefore there are at least two rows of length 
> e 2 N in the Young tableaux, and hence the dimension of the corresponding 
representation is at least the one of the tableaux consisting of two rows of 
length I = [e 2 N] (here the square bracket denotes the integer part). This is 
given by 

(2f)l &*■ 

WTTy-W («2) 

which proves the Corollary. 

To prove the statement on the upper bound on the ground state degener- 
acy in Remark 3 after Cor.QJwe estimate, for 6(z) = 1 for z > 1 and 9(z) = 
for z < 1, 

jV(N, \nigN~ 1 ) = TrleiE^iN^gN-^-H^jf-i)] 

< Tr[6(N(E DM (n,g) + e)-ZlX l) )] (4-83) 

for any e > and N large enough. Here we have used ()3.3|) and the fact that 
the interaction potential v is positive by assumption. Estimating the step 
function by the exponential function, we therefore get, for any j3 > 0, 

Af^N^nigN- 1 ) < Tr[evp((1(N(E™(n,g)+e)-j:l 1 H®))] 

= Tr[exp{(3(E DM (n,g)+e-H ))} N , (4.84) 



which proves (j3.13|) . 
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A Appendix: Proof of Theorem U 

The strategy of the proof of Theorem ^ is similar to the one in the two- 
dimensional case given in in [TT]. One shows the non- uniqueness of mini- 
mizers of the GP functional by showing that no eigenfunction of the angular 
momentum f2 • L can be a minimizer, even though the GP functional is in- 
variant under rotations around the f2 axis. This implies that the minimizer 
cannot be unique, since by rotating a minimizer one gets another one. More- 
over, since the minimizer of the DM functional is unique by Prop. Q 
this implies that it has to have at least rank 2, and E DM (Q,g) < E GP (tt, g), 
as stated in Theorem ^ 

Thus all we have to show is that an eigenfunction of f2 • L cannot minimize 
£fi g- To do this, we first have to consider the properties of minimizers in the 
subspaces of angular momentum eigenfunctions. Except for Lemma HI all 
the considerations below follow the same line as in the two-dimensional case. 
Lemma IHl however, is the main new ingredient to extend the results to three 
dimensions. While Lemmas OEI hold for general V's, our assumptions (|2.8|) 
and ()2.9|) on the external potential V will become important in the proof of 
Theorem [T] below, where they are used in order to get explicit estimates on 
the various quantities appearing in Lemmas EHZI 

Let f n denote the minimizer of the GP functional restricted to functions 
<f)(x) = f(r,z)e mip , i.e., to functions with angular momentum f2 • L = n. 
More precisely, f n minimizes 



f~ d 2 r dzf* -A r f --L+ f + V f + 47rg\f\ 2 f (A.l) 



under the normalization condition J \f(r, z)\ 2 d 2 rdz = 1, where we denoted 
d 2 r = 2irrdr, and A r = d 2 /dr 2 + r~~ 1 d/dr. In the following, it will be 
convenient to study this functional for all n e M + , not only for integers. 
Using standard methods one can show that for each n > and g > there 
exists a unique minimizer of (|A.1|) . Let E n (g) denote the corresponding 
minimum of (jA.l|) . and let 






(A.2) 



Note that jj, n depends on g besides n. 

We need the following estimates on the minimizers /, 
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LEMMA 5. The minimizer f n is a non-negative, bounded function, with 

ll/»IU < ^ • (A.3) 

Proof. As in Lemma 1] this follows from the maximum principle, applied 
to the variational equation for /„, which reads 

-AJ n - —f n + —f n + Vf n = (Jl n - 87rg\f n \ 2 ) f n . (A.4) 

Recall that V > by assumption. □ 

LEMMA 6. Let j n (r) = j°° oo dz\f n (r ) z)\ 2 . Choose R > large enough such 
that 

5 n (R) = inf V(x) -Jln>0, (A.5) 

\z\>R 

which is always possible by our assumption that V goes to infinity at infinity. 
Then 

\\jn\\oo< || /n|| l(^R+ COnst. • (A.6) 

Here \\ ■ \\ has to be understood as the supremum norm on M. 2 and M 3 , re- 
spectively. 

Proof. Let r] e C°°(R) be such that < rj < 1, rj{z) = 1 for \z\ > 2 and 
r)(z) = for \z\ < 1. Then 



Jn (r) = / dzlU^z^dl-^z/R^+^z/R)) <4R\\f n \\l + X (r) , (A.7) 
where we denoted 

X (r)= [ dz\f n (r,z)\ 2 V (z/R) . (A.8) 



Using the variational equation for f n and neglecting positive terms, it is 
straightforward to derive the inequality 

-A rX < 2 / dzrj(z/R) (\f n \ 2 (Ji n - V(r, z)) + f n -^fn) ■ (A.9) 
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Partial integration yields the estimate 

2 J dz V (z/R) (/n^/n) <^J dz\f n \Y( Z /R) < C -^\\f\\l = S . 

(A. 10) 

Using (jA.5|) and the support properties of r\ we therefore get 

-A rX <e-25 n (R) X . (A.ll) 

It follows from the maximum principle that x( r ) ^ e l (2(5„(i?)) , which proves 
our claim. □ 

LEMMA 7. With j n given in Lemma\^ r~ 2n j n (r) is a bounded function, 
with 

||r- 2n j„||oc < (2c^„) n ||j„,|| 00 , (A.12) 

where c n is given by 



' 2 -«(2_«)" /2 7 rCsc(f) xl/n 



for n < 1 



n Tin 



(2 - n)r(n) 

forn>\. (A. 13) 



Proof. The proof of this Lemma is similar to the corresponding Lemma 2 in 
|llj . Let /i(r) = r~ 2n x{r). Using the variational equation for f n and the fact 
that V is positive by assumption, one derives the differential inequality 

-h"( r ) - ^ 2n+l \ '(r) < 2a n h(r) . (A.14) 
r 

Multiplying (jA.14|) with the integral kernel of the resolvent of the operator 
on the left side of this inequality and integrating, we can proceed as in [TT| 
Lemma 2] to arrive at (jA.12|) . □ 

We now are able to give the proof of Theorem ^ 

Proof of Theorem^ We start by estimating the difference of E n+1 (g) and 
E n (g). By analogous considerations as in JTTJ Eqs. (3.2) and (3.7)-(3.11)] we 
have, using Lemma [7| in replacement of Lemma 2], 



E n+1 {g) - E n {g) < (2n + 1) [E^g) - E (g)) 

+ ln ( 

0<m<l L 



<(2n + l)2vre sup ||j m ||oo max \ 3, 1 + In ( — f m ) \ . (A.15) 
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We now use the bounds (|A.6J) and (|A.3|) in order to estimate ||j m ||oo- Our 
assumptions on the external potential V imply that, if S n (R) = 0(1), 

~ 9 1/{1+2/s+1/p) and R~^ n /P (A. 16) 

for large g, uniformly in n for bounded n. Therefore 

E n+1 (g) - E n {g) < O L^W*) Info)) ^ . (A.17) 

Note that the minimum of the GP functional over all functions with angular 
momentum n is given by E n (g) — n\fl\, and hence, for any fixed given N, all 
functions with angular momentum < A" cannot be GP minimizers if 

min {E n (g) - n\fl\} > min{E n (g) - n\fl\} . (A.18) 

0<n<N n>N 

This is fulfilled if 

IOI ^ En ^) ~ E n{g) /a -• q\ 

S2 > max , (A. 19) 

0<n<AT N — n 

which holds true, for any given fixed A^ and f2 7^ 0, for large enough g by 
(|A~T7l) . 

It remains to show that for given f2 there exists an A^, independent 
of g, such that all functions with angular momentum > cannot be GP 
minimizers. We do this by showing that they are unstable, in the sense 
that a small perturbation lowers the energy. This can be done by a similar 
calculation as in the two-dimensional case [TI], using the estimates of Lem- 
mas |3HZI above. We will only sketch the argument here, mainly pointing out 
the difference to the two-dimensional case. 

In order to show instability, one has to show that the second derivative 
of the GP functional at some given stationary state 4>{x) = f n (r,z)e vnip is 
negative in some direction. More precisely, 

Q(w) = (w\H + l6irg\f n \ 2 -Jl n + n\n\\w) + 8ng$l J (j)* 2 w 2 < (A.20) 

for some w with (<f)\w) = 0. Here 5i denotes the real part of a complex 
number. We will choose an w that, when added to the function f n (r, z)e tnip , 
effectively splits the central vortex of degree n into n different vortices of 
degree 1 that are arranged symmetrically around the axis r = 0. The idea is 
that for large enough n a splitting of a central vortex like this should lower 
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the energy of the GP functional, and therefore the second derivative given 
by Q(w) should be negative. 

Let n > 1, and let W\ G H l (M?) be a radial function, with || 1 1| 2 = 1 and 
support in a ball of radius 1. Let v G iJ 1 (M), with ||t>|| 2 = 1, and choose as 
a trial function for Q 

w(r,z) = (2c^ n )VV((2c^„) 1/2 r) I ^,;(z/L) , (A.21) 

with c n given in (jA.13|) . and L > 1. Note that c n = 0(n~ l l 2 ) for large n and, 
by (HH and (lA~T|) . 

Jin > inf + V(x)\ > noo - C w (A.22) 

a; |^ r 2 J 

for any uo < Q c and some constant depending only on u. Therefore 
c^fin > const, independent of n and g. Using ()2.8j) this implies that 

(w\V\w) < const. L p . (A.23) 

Denoting T r = (wi\ — A r \wx), T z = (v\ — d 2 /dz 2 \v) and M n = j d 2 rr 2n wi, 
we get the bound (note that the last term in ()A.20)1 is zero for the w under 
consideration, since n > 1) 

Q(w) <n|0|-/I n +2c^ n T r +T,+const.L p +167r|||r- 2n j n (r)|| 00 M n . (A.24) 

For some fixed e > and n large enough, choose L such that the right side 
of (jA.23|) equals ejl n . For this choice of L, 

\\r~ 2n j n \\oo < C £ -Jl n (A.25) 
9 

for some constant C e depending only on e, as can be seen from Lemmas EH3 
and the fact that, for large Jl n , R ~ Jm! p if S n (R) = 0(1) in (jA.5|) . again by 
our assumptions ()2.8|) and ()2.9|) on V. Using these estimates in (jA.24|) we 
obtain 

Q(w) < n\n\ -Ji n {l-e- 2c 2 n T r - 16ttC ; £ M„,) + T z . (A.26) 

Since |f2| < Q c by assumption, and both (? n and M n tend to zero as n — » 00, 
we can choose e appropriately to conclude from (|A.22|) that the right side of 
(jA.26|) is negative for n large enough. 
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We have thus shown that for some Nfi, independent of g, all minimizers 
of the GP functional restricted to eigenfunctions of the angular momentum 
f2 ■ L with eigenvalue n > are unstable and therefore cannot be absolute 
minimizers of the GP functional. Together with the result (jA.19|) above this 
proves Theorem □ 
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